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Abstract 

The generation of ultra-strong attosecond pulses through laser-plasma interactions offers the 
opportunity to surpass the intensity of any known laboratory radiation source, giving rise to new 
experimental possibilities, such as quantum electrodynamical tests and matter probing at extremely 
short scales. Here we demonstrate that a laser irradiated plasma surface can act as an efficient 
converter from the femto- to the attosecond range, giving a dramatic rise in pulse intensity. Al- 
though seemingly similar schemes have been presented in the literature, the present setup deviates 
significantly from previous attempts. We present a new model describing the nonlinear process of 
relativistic laser-plasma interaction. This model, which is applicable to a multitude of phenomena, 
is shown to be in excellent agreement with particle-in-cell simulations. We provide, through our 
model, the necessary details for an experiment to be performed. The possibility to reach intensities 
above 10 26 W/cm 2 , using upcoming 10 petawatt laser sources, is demonstrated. 

PACS numbers: 52.38.-r, 52.27.Ny, 42.65.Ky, 42.65.Re 
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Introduction. Recent progress in u 
pulse intensities surpassing 10 22 W/cm 2 



trahigh-power laser technology has resulted in 
and stimulated the construction of multi- 



petawatt laser sources 0. Sueh .asers open up opportuuit.es for studying both a number 
of fundamentally new problems, such as the effects of vacuum nonlinearities [3H5( in laser 
fields and photo- nuclear physics, as well as some very important applications, e.g., laser 
based particle acceleration, fast ignition fusion schemes, and the generation of electromag- 
netic radiation with tailored properties. 

Given this, the study of overdense plasmas irradiated by relativistically intense laser 
pulses constitutes a very important and challenging research direction. Numerical studies 
using particle-in-cell approach, taking into account most important effects for the typical 
range of parameters, are known to be an excellent tool in this field, and the numerical results 
in general agree well with the experimental results. Furthermore, so called nonlinear fluid 
models jg] give a set of equations analytically describing such processes, but the strongly 
nonlinear plasma behavior, due to the ultrarelativistic motion of the plasma electrons, makes 
the development of a theoretical approaches a highly complex task. Thus, one is normally 
forced to limit oneself to qualitative analyses, and use a phenomenologically motivated ad 
hoc treatment. 

The generation of high harmonics from intense laser-plasma interactions is an intensely 
studied research field, with a manifold of applications pfl, including the idea of reaching the 
extreme intensities needed to probe vacuum nonlinearities using lasers j^- 12|. 

As of today, the most prominent theoretical model used in the analysis of such high- 
order harmonics generation (HHG) is the so-called oscillating mirror model (OMM). In the 
OMM, one considers the backradiation from an overdense plasma by taking into account 
the retarded emission from the oscillating source. This approach was first proposed by 

Il5l |. Lately, the OMM approach 



Bulanov et al. [131 and developed further in Refs. 14L 

n 

has been reexamined by Gordienko et al. [lq . who proposed that at each moment of time 
there exists a so-called apparent reflection point (ARP) at which the energy flux vanishes. 
This assumption implies a local (in time) energy conservation or, phrased differently, the 
approach neglects the energy accumulated by the plasma, in the form of the fields due 
to charge separation caused by the light pressure. An asymptotic analysis of the ARP 



dynamics in the strongly relativistic limit 17j indicates the universal properties of the HHG 
spectra: the intensity of nth harmonic scales as n -8 / 3 , and the cutoff ~ 7 max , where 7 max is 
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FIG. 1: The parameter 5 obtained from ID PIC simulations of a plasma with density N e at oblique 
irradiation (0 = 60°, P-polarization) by a wave with constant intensity I and a 1/im wavelength 
during one optical cycle. {S is the ultrarelativistic similarity parameter, defined as the quotient 
between dimensionless density and the dimensionless amplitude (see Eq. (2)).) 

;he maximal relativistic factor of the ARP. These results agree well with the experiments 
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20]. 



Nevertheless, the assumption of the OMM concerning local temporal energy conservation 
is valid only for restricted values of the plasma density and laser intensity. The parameter 
5, used in Fig. [TJ is defined as the ratio between the maximal accumulated plasma field 
energy and the energy of one optical cycle. In the bottom left corner of the figure we find 
the zone labeled "OMM", for which 5 <C 1. Thus, here the energy accumulation can be 
neglected, as assumed in the OMM. In the top right corner the region of relativistically 
self-induced transparency (RSIT) is shown. Thus, there is a large, and very important, 
parameter region that so far has not been covered by any theoretical model, and for which 
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5 ~ 1 such that the OMM's assumption of local energy conservation is no longer valid. Here 
collective ultrarelativistic electrons motion can give rise to a nanoplasmonic structures, i.e. 
nanometer scale surface layers and their oscillations, absorbing at each period the energy 
of the laser pulse to the energy of internal electric and magnetic plasma fields and then 
reemmiting this energy in the form of attosecond burst. The process result in a markedly 
slower decay in the generated higher harmonic spectra 2l| as compared to the OMM results. 
Furthermore, based on the phenomenological assumption of electron nanobunches appearing 
in the plasma, emitting radiation, it was shown in Ref. 22j that the spectra can be much 
flatter than predicted by the OMM. Moreover, the OMM assumes as a prerequisit that the 
incident and backradiated amplitudes at the ARP are equal, in accordance to the Leontovich 
boundary conditions (which is in direct correspondence to the local energy conservation). 
Consequently, situations where large field amplifications is to be expected cannot be analyzed 
through this model. Thus, finding a new theoretical model in the relevant parameter regime, 
between the OMM and RSIT regions, is of the utmost importance for a large number of 
applications. 

In the present work we propose, for the first time, a physically motivated model, a so- 
called relativistic electronic spring (RES) model, describing the highly nonlinear behavior 
of laser-plasma interactions. The model gives very good agreement with simulation and 
makes it possible to analytically study a vast range of regimes in laser-plasma interactions 
that otherwise would be out of reach for analysis. In particular, one of the most remarkable 
effects in the RES regime is the possibility of generating attosecond pulses with an amplitude 
several orders of magnitude higher than the incident laser pulse. Here, we apply the RES 
model to this amplification effect in order to understand the underlying physical mechanisms 
and determine the optimal parameters for an experiment to be performed. We compare the 
results with particle-in-cell simulations and find excellent agreement. The implications of 
our results are discussed, in particular the possibility to utilize this new type of secondary 
source for novel experiments. 

Ultrarelativistic energy conversion on the surface of a planar target. - The process of 
giant attosecond pulse generation due to the oblique incidence of a p-polarized electromag- 
netic wave on a plane plasma boundary (x = 0, y, z) may be considered using a boosted 
frame moving along the plasma surface and plane of incidence in the y-direction with the 
velocity csin6*, where c is the speed of light and 6 is the angle of incidence (see Fig. [21(a)), 
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thus making the problem one-dimensional 23]. 

The incident laser pulse pushes the electrons into the plasma due to the light pressure. 
Unlike the case of normal incidence, the oblique incidence results in the emergence of un- 
compensated currents and magnetic fields in the boosted frame. Therefore, the electrons 
experience an additional ponderomotive action, which is different during the two half-periods 
of the incident wave. When the laser electric field is directed along the y-axis, the Lorentz 
force due to uncompensated currents enhances the light pressure effect, pushing electrons 
further from the boundary. It can be analytically shown that under the laser radiation 
pressure the electrons are shifted and group into a boundary layer. The thickness of this 
layer tends to zero as the laser intensity increases and, in the ultrarelativistic limit, is much 
smaller than all other spatial scales involved in the process. The charge and current densities 
in this layer greatly exceed the ones in the unperturbed plasma. At this stage the incident 
wave energy is transformed to the energy of the internal plasma fields and kinetic energy of 
the particles. 

The formation of an ultrathin (nanoscale) electron layer due to the interaction between 
an ultrarelativistic laser pulse and an overdense plasma has been known for about a decade, 
being reported in works on relativistic self-induced transparency 24], |25] and particle accel- 



eration using thin foils [26]. Unlike in the case of a circularly polarized laser pulse, in which 
electrons may be shifted by a distance of several wavelengths, a linearly p-polarized pulse 
results in a light pressure force that oscillates within a field period. Hence, the electrons 
are pushed from the surface for no longer than a fraction of the optical wavelength, then 
break away under the action of the charge separation force and travel towards the incident 
wave in the form of a nanobunch, providing a source of attosecond burst. The presently 
described processes can clearly be seen in Fig. Mjo), where the results of ID PIC simulation 
are presented. 

Thus, the described process and concomitant energy conversion may be represented as a 
sequence of three stages: 1) the pushing of electrons from the surface by the ponderomotive 
force and the formation of a thin current layer giving an energy transfer from the laser field to 
the plasma fields and particles; 2) the backward accelerated motion of the electrons towards 
the incident wave with the conversion of the energy accumulated in the plasma and laser 
field energy into kinetic energy of an ultrarelativistic electron bunch; and 3) the radiation 
of attosecond pulses by an electron bunch due to conversion of the kinetic energy and laser 
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FIG. 2: (a) Transformation to a moving reference frame. Here, the plasma density n, wave ampli- 
tude a, and frequency to have a subscript denoting the moving frame, (b) Space-time distribution 
of electrons density N (green) , amplitudes of incident cij n (red) and backradiated a g (blue) electro- 
magnetic fluxes obtained from ID PIC simulation of plasma with density 4 x 10 23 cm -3 at oblique 
irradiation (9 = 11.25°) by a wave with constant intensity 10 23 W/cm 2 and 1 fim wavelength 
during three optical cycles. The red curve shows position of the thin layer obtained using the RES 
model. Time and coordinate are in dimensionless units (see. Eq. ([2])). (c) Backradiated signal 
obtained by PIC simulation (blue) and using the RES model (red). 



field energy to the XUV and X-ray range. Based on the motion of the plasma electrons and 
the energy conversion scenario we find that it is natural to refer to this three-step process 
as to a model of a relativistic electronic spring (RES). It should be emphasized that due 
to the energy accumulation in the plasma, the backradiated field can be much larger than 
the incident field. This is the fundamental difference from the OMM (see, e.g., 27, 2sj|), a 
model which assumes a direct correspondence between the incident and backradiated field. 
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The RES model. A model describing the dynamics of a thin current layer and the 
generation of attosecond pulses may be formulated starting from three intuitively clear and 
physically justified prerequisities, that can be verified using PIC simulation. First, we will 
assume that at each moment of time the plasma electrons are represented by two fractions: 
one infinitely narrow layer of shifted electrons at a certain moving point x s , where all the 
electrons from the region < x < x s are accumulated, and one population of electrons with 
unperturbed density at x > x s . Second, we will characterize the motion of particles in the 
layer by the ensemble averaged x- and y-projections of the velocities f3 x and (3 y (normalized 
by the speed of light, and with the 7-factor defined by 0%. + f3y = 1 — 7 -2 ). Third, we 
will suppose that the motion of the electrons in the thin layer together with the flow of 
uncompensated ions in the < x < x s region fully compensate the incident electromagnetic 
radiation in the unperturbed plasma at x > x s . 

It is readily shown from Maxwell's equations that in a one-dimensional geometry, a 
moving charged layer, with surface charge a, emits electromagnetic waves with amplitudes 
2naf3 y / (1 — /3 X ) and 2^af3 y / (1 + j3 x ) in the positive and negative directions of the rr-axis, 
respectively. Consequently, the expression for the incident wave compensation motion may 
be written 

where the left-hand side corresponds to the incident wave, whereas the terms on the right- 
hand side describe the radiation of the uncompensated ions and the electron layer, respec- 
tively; here S = n/a is the ultrarelativistic similarity parameter, while a and n are the 
incident wave amplitude and the plasma density in the laboratory frame, respectively. We 
use dimensionless quantities which can be expressed in terms of dimensional time t, coordi- 
nate x, density n and electric field amplitude a according to 

~ tocos 9 „ Aire 2 A e A 

t = u cos 9t, x = x, n = -n, a = a, (2) 

c mar mau 

where u is the carrier laser frequency, while m and e are the electron mass and charge, 
respectively. Analogous to the Eq. (JTJ, the electric field, as a function of retarded time 
£ = x + t, emitted by the plasma in the negative x-direction is given by 

a 9 g = x s (t) + t] = a 0J ^- e - sin9) z.(t), (3) 

where a = a cos 9 is the incident wave amplitude in the boosted frame. 



The layer dynamics is determined by the equation 

j t x a = P. (4) 

with the initial condition x s (t = 0) = 0. By virtue of the ultrarelativistic motion, the 
position of the layer may be found by assuming that the full particle velocity is equal to the 
speed of light, i.e., /3 2 + /3 2 = 1. Equations (CD) and (HJ are then self-consistent and the layer 
motion is described by a first-order nonautonomous ordinary differential equation or by an 
autonomous system 

2 _ i N Q _ „.\ _|_ A ^o3 Q C-l h _ „2 _ „,\2H X / 2 



rf n (M 2 -l)(sin^-M)±4cos 3 ^5- 1 l- 77 2 (sin# 



rfr 7] (u 2 + 1) (5) 

drj u 2 — 1 

. dr u 2 + 1 

for the variables 77(f) = x s 5y (2cos 3 #), u(r) = /3 y / (1 — (3 X ), where r = tS/ (2cos 3 8). The 
solution of Eq. (JH]) depends on the two dimensionless variables 5 and 9, and may be analyzed 
in the plane {77, u}, where we have two sheets corresponding to the choice of sign in Eq. flHJ). 

The topology of the phase plane is characterized by the existence of a stable limit cycle 
(see Fig. [3]). It is convenient to classify the form of the solution by the number of zero-axis 
intersections of u{r) or, equivalently, by (3 y (t) in the 77 > (x s > 0) region. This takes 
place when (3 X — > — 1 and Eq. ([3]) becomes singular, which corresponds to the emission 
of the attosecond burst. There may occur either two such events (Fig. GJc)), one (Fig. 

or none (Fig. [3(d)) in each optical period. Accordingly, two bipolar, one bipolar, or 
one unipolar attosecond pulse are generated. It is clear that we have the emission of two 
evenly spaced identical bursts in the case of normal incidence (Fig. 02(b)). With increasing 
angle 9 the second burst either disappears due to amplitude decay down to zero ((c) —> (e) 
transition) or two bipolar pulses merge into one unipolar pulse as a result of convergence of 
their generation times ((c) — » (d) bifurcation). Our comprehensive numerical study indicates 
that the results obtained using the RES model are in a very good agreement with the PIC 
simulations for all values of 9 and for S < 5 down to RSIT (in particular, see Fig. [2(b), (c)). 
The RES model's applicability region is labeled "RES" in Fig. [TJ 

Note that by a simple modification of Eqs. ([T|) and ([3]), the RES model can be easily 
generalized to take into account an arbitrary plasma density profile, as well as arbitrary 
laser pulse shape and polarization. 
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FIG. 3: (a) Regions on the plane of parameters S, 9 corresponding to qualitatively different forms of 
solution of Eq. ([5]); (b) — (f) form of solution on the plane {rj, u} (red), its limit cycle (black), lines 
of phase portrait for the + sign (blue) and the — sign (green) in Eq. ©, and the corresponding 
form of electromagnetic backradiation from the plasma. 

In order to find the amplitude and duration of the pulse generated near j3 y = one has to 
take into account the finite value of the relativistic factor 7, which is the external parameter 
for the RES model and can be taken, for example, from PIC simulation. Near the point 
(3 y = 0, /3 X — > — 1 the first term in the right-hand side of Eq. ([3]) dominates, so an analytical 
expression for the burst shape and its spectrum can be written as follows 



where a = ^j- 



A,=cr 



S(0 = A 9 f 



h oc exp 



cry- 



TO' A 9 = a ° S 2co% V ° is the P ulse amplitude, r g 



[2a~f 



(6) 

3V-1 



is its characteristic duration, and 1^ is the intensity of the kth harmonic. The value of 
a is assessed from the solution of the self-consistent system; hence, it depends only on 
dimensionless parameters S and 9 and is independent of 7. This means that pulse duration 
in the ultrarelativistic limit tends to zero as 7 -3 . Note that in the RES model the spectrum 



9 



decays exponentially with the characteristic scale cry 3 and, in contrast to the OMM, has no 
region with a power-law decay. Thus, the "RES" region in Fig. [1] corresponds to a slower 
energy decay of the harmonics than the "OMM" region. 

Depending on the thickness l s of radiating electron layer, the radiation can be either 
coherent, for l s < r g , or incoherent, for l s > r g . In the latter case, the assessment of the giant 
pulse amplitude requires taking into account that only part of electrons radiates coherently. 
The radiating electron layer thickness l s may be estimated assuming that, starting from 
the time of maximum displacement, it decreases proportionally to the number of electrons 
in the layer. Besides, as the electron velocity relative to the £ coordinate decreases from 
the speed of light to c(l + /3 X ), the layer experiences an additional 27 2 -fold compression. 
The layer thickness L s at maximum displacement may be estimated from the equation of 
balance between the force of light pressure and forces caused by charge separation in the 
form L s fa y2 cos 9S~ 1 ^ 3 a~ 2 ^ 3 . The resulting estimation of the adjustment factor is 

c = r f = x4— x r* V 1/3q2/3 ' ( ? ) 

l s v 2 «7 x s\p y=0 cos # 
where x max is the maximum value of x s , and \ is a dimensionless constant of order unity that 
is needed to account for arbitrary choice of the estimated values entering this expression. 

The PIC simulations and the results of the numerical solution of Eqs. ([!]), 01]), including 
the adjustment factor ([7]), are compared in Fig. HI We see that we have not only qualita- 
tive, but also quantitative agreement for the burst amplitude and the generation phase. A 
qualitative agreement can be seen for the burst duration and regions of unipolar generation. 
The diagrams in Fig. HJ^b, c) allow us to distinguish the zone with the center 

fl «62°, S g ^ 1/2, (8) 

and the boundaries 50° < 9 < 70°, l/4<S'<las the region of the most powerful 
and short burst generation. This result may serve as a guiding message for experimental 
implementation. The optimal parameters (jHJ) correspond to the triple point in the S and 
9 plane (Fig. [3^a, f)), which provides the longest time of irradiation and thus maximum 
amplitude of the generated attosecond pulse, as coherency plays the dominant role. 

A focusing mechanism for attosecond pulses. - Based on the results obtained we pro- 
pose a new concept of extremely intense light generation at the level required for observation 
of the QED effects. The idea is to focus the giant burst formed in the regime described above 
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FIG. 4: Generation phase £ s , maximum amplitude a g /ao and duration r g of giant pulse obtained 
from the RES model with 7 = 30 and x = 1/5 (d, e, f) and from ID PIC simulations (a, b, c) 
of semi-infinite plasma obliquely irradiated by one optical cycle pulse with 1 /j,m wavelength and 
amplitude a = 191.1, which corresponds to 10 23 W/cm 2 intensity. For the PIC simulations, the 
pulse duration was assessed as a distance between maximum and minimum electric field points, and 
the region of unipolar pulse generation is given in white. On the (f) panel, the coherent generation 
and monopolar pulse generation regions are marked (1) and (2), respectively. 

by using a slightly grooved surface of the obliquely irradiated target at the optimal parame- 
ters (JHJ), with the guiding line of the groove located in the plane of incidence (see Fig. [3(a)). 
The PIC simulation of the proposed concept shows that the intensity 1.8xl0 26 W/cm 2 can 
be reached in the zone with the size of order 10 nm with a 10 PW laser pulse, as can be seen 
in Fig. UJb). In the laboratory frame the high field zone moves along the guiding line with 
the speed c/ sin#. We note that the size of laser pulse along the transverse direction may 
be only a few wavelengths. The data shown in Fig. H]may be used to modify surface profile 
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FIG. 5: (a) Schematic representation of the concept of groove-shaped target, (b) Intensity dis- 
tribution at focusing instant obtained from 2D PIC simulation: a linearly polarized wave with 
intensity 10 23 W/ cm and wavelength 1 /im obliquely incident at an optimal angle 6 g = 62° on a 
parabolic groove-shaped target with density 1.065xl0 23 cm" 3 , which corresponds to S = 1/2. 

and target density to allow using a laser pulse with a more complicated intensity profile in 
the transverse direction. 

Conclusion. In this work we studied the giant pulse generation process at oblique 
irradiation in an overdense plasma by a relativistically strong laser pulse. The model of 
relativistic electronic spring was developed, providing a qualitative and, for some charac- 
teristics, also a fairly good quantitative description. The parameters of the most powerful 
burst generation (jSJ) were determined. A new concept of a groove-shaped target for high 
electromagnetic field generation aimed at obtaining the QED effects by means of upcoming 
laser sources was proposed and confirmed by PIC simulation. 

Methods. For the numerical studies we used a particle-in-cell approach. The data 
shown in Fig. (OQ), (j2J) and (J4]) were obtained using ID simulations in the moving frame (to 
take into account the oblique incidence of the laser pulse), implying a plasma flux in the 
transverse direction. For Fig. (j3j) we have carried out a 2D simulation in the frame moving 
along the guiding line of the groove. Performing the PIC simulation in the moving frame 
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implies that the transverse size of the laser pulse is fairly large compared to the wavelength: 
However, this is a rather weak restriction. A fully relativistic parallel FFT based PIC code 
("ELMIS" 29] ) have been used; an 8 /im x 8 /xm region is represented by 8192x8192 cells. 



The plasma ions are taken to be Au 6+ , and each target cell contains approximately 100 
virtual particles of each species. The time step is 1.07xl0~ 2 fs, and the laser pulse front has 
a sine-squared profile with a two wave periods duration. 
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